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(57) Abstract 

In order to cancel echo on several sound signal channels, the 
method consists in generating on each channel a synthetic signal having 
the spectral characteristics of a white noise; computing the frequency 
masking thresholds in several adjacent frequency bands using auditory 
perception properties; bringing in each frequency channel the synthetic 
signal level to the associated masking threshold value, so as to obtain 
an auxiliary signal, which is added to the signal associated with the 
loudspeaker of the channel concerned; supplying the resulting signal to 
an adaptive filter, adding the output signals of the associated adaptive 
filters to each microphone channel; subtracting the resulting signal 
received by the microphone of this channel; deducing an estimation 
error, on the basis of which the coefficients of the adaptive filters 
associated with this channel are corrected. 

(57) Abrege* 

Pour effectuer une annulation d'echo sur plusieurs voies de signal 
sonore, sur chaque voie, on engendre un signal synth6tique ayant les 
caracteristiques spectrales d'un bruit blanc; on calcule des seuils de 
masquage frequentiel dans plusieurs bandes de frequence adjacentes, 
en utilisant des proprtetes de la perception auditive; dans chaque bande 
de frequence, on amene le niveau du signal synth£tique a la valeur du 
seuil de masquage associe\ de facon a obtenir un signal auxiliaire, qu'on 
additionne au signal associe au haut-parleur de la voie consid6r6e; on 
fournit le signal obtenu a un filtre adaptatif; on additionne les signaux 
de sortie des filtres adaptatifs associes a chaque voie de microphone; 
on soustrait le signal resultant du signal recu par le microphone de cette 
voie; on en dSduit une erreur d'estimation, en fonction de laquelle on 
corrige les coefficients des filtres adaptatifs assocfes a cette voie. 
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PROCEDE D f ANNUL ATI ON D ' ECHO ACOUSTIQUE MULTI-VOIES ET 
ANNULEUR D T ECHO ACOUSTIQUE MULTI-VOIES 

La presente invention concerne un procede d ' annula- 
5 tion d'echo acoustique multi-voies, ainsi qu'un annuleur 
d'echo acoustique multi-voies. 

Dans le domaine de la transmission des signaux 
sonores, dans certaines applications telles que la telepho- 
nie "mains libres" et la teleconference, l'echo acoustique 
10 est une source de gene considerable. Les dispositifs connus 
qui combattent 1 1 echo acoustique comportent generalement des 
f litres adaptatif s, dont le role est d ? identifier et de 
modeliser la reponse impulsionnelle du chemin de couplage 
acoustique entre le ou les haut-parleurs et le ou les 
15 microphones du systeme de transmission de signaux sonores 

considere. 

La figure 1 rappelle la structure generale d'un 
annuleur d'echo acoustique classique. II est associe a un 
haut-parleur 10 et un microphone 12 entre lesquels il existe 

20 un chemin de couplage acoustique ou chemin d'6cho 14. Les 
references 16 et 18 designent respectivement le signal 
sonore regu et le signal sonore emis. L' annuleur d'echo de 
la figure 1 comprend un filtre adaptatif 20 recevant en 
entr6e le signal sonore re?u 16. L'objet du filtre adaptatif 

25 20 est d 1 estimer par ses coefficients la reponse impulsion- 
nelle du chemin d'echo 14, afin de soustraire l'echo du 
signal regu par le microphone 12. Pour cela, la sortie du 
filtre adaptatif 20 est connectee a un soustracteur 22 qui 
soustrait le signal obtenu en sortie du filtre adaptatif 20 

30 du signal regu par le microphone 12. Le signal de difference 
obtenu en sortie du soustracteur 22 fournit une erreur 
d' estimation 24. Les coefficients du filtre adaptatif 20 
sont ajustes au cours du temps par un algorithme approprie 
qui exploite 1 ' information d' erreur d ' estimation. 

35 Le choix de cet algorithme conditionne de fagon 

determinante les performances de 1 ' annuleur d'6cho. Les 
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procedes et dispositifs d 1 annulation d'echo connus utilisent 
couramment un algorithme dit du gradient stochastique 
normalise, habituellement designe par le sigle NLMS. Un 
inconvenient: de cet algorithme est que sa vitesse de 
5 convergence depend des caracteristiques spectrales du signal 

sonore regu. Un objectif de la presente invention est de 
reduire cette dependance. 

Le probleme de l'echo acoustique se pose bien sur 
egalement dans les systemes multi-voies de transmission de 
10 signaux sonores, c'est-a-dire dans les systemes a plusieurs 
haut-parleurs et plusieurs microphones, par exemple deux 
haut-parleurs et deux microphones dans le cas de la stereo- 
phonie. 

On connait des procedes et des dispositifs d'annula- 

15 tion d'echo multi-voies fondes sur le meme principe que ceux 

d' annulation mono-voie. La figure 2 represente a titre 
d' exemple le schema de principe partiel d'un annuleur d'echo 
acoustique stereophonique classique, ou, pour plus de 
clarte, on n'a represente qu'une seule des deux voies de 

20 microphone. La structure generale peut etre ai semen t 

generalisee a celle d'un annuleur d'echo acoustique a N 
voies de signal sonore, ou N est un entier quelconque. On ne 
decrit dans ce qui suit que le cas N = 2. De fagon analogue 
a 1' annuleur d'echo de la figure 1, 1 ' annuleur d'echo 

25 stereophonique s' applique a des voies de signal sonore regu 

16 x , 16 2 et des voies de signal sonore emis, dont une seule, 
designee par le chiffre de reference 18^ est representee. 
L' annuleur d'echo est associe a deux haut-parleurs 10 a/ 10 2 
et deux microphones, dont un seul, designe par le chiffre de 

30 reference 12 x , est represente. II existe quatre chemins 
d'echo, dont deux chemins d'echo (14 x et 14 2 ) entre les deux 
haut-parleurs 10^ 10 2 et le premier microphone 12 1 et deux 
autres chemins d 1 echo ( non representes ) entre les deux haut- 
parleurs 10 lf 10 2 et le second microphone. Pour estimer les 

35 reponses impulsionnelles des divers chemins d'echo, un 
f litre adaptatif est prevu entre chaque voie de haut-parleur 
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et chaque voie de microphone. Ainsi, si on considere la voie 
du microphone 12 x , il est prevu deux filtres adaptatifs 20 lf 
20 2 qui regoivent respectivement en entrfee les signaux 
sonores regus 16 x et 16 2 . Les sorties des filtres adaptatifs 
5 20 lf 20 2 sont fournies en entree a un additionneur 26 t . Dans 
un annuleur d'echo a N voies.ou N est un entier quelconque, 
on prevoit un tel additionneur sur chaque voie de micropho- 
ne. La somme des signaux de sortie des filtres adaptatifs 
fournie par 1 ' additionneur 26 x est soustraite, par un 
10 soustracteur 22 x , du signal regu par le microphone 12 x . La 
meme operation est effectuee sur chaque voie de microphone. 
Le signal de difference obtenu en sortie du soustracteur 22 x 
fournit une erreur d' estimation 24 x commune a tous les 
filtres adaptatifs relatif s a la voie de microphone conside- 
rs ree ; dans l'exemple de la figure 2, l 1 erreur d 1 estimation 
24 x est commune aux filtres adaptatifs 20 x et 20 2 . De meme 
qu'en annulation d'echo mono-voies, les coefficients des 
filtres adaptatifs sont modifies de fagon iterative par un 
algorithme approprie, sur chaque voie de microphone, & 
20 partir de 1 1 erreur d f estimation obtenue. 

On a constate, en annulation d'echo multi-voies, et 
en particulier en annulation d'echo stereophonique, que la 
vitesse de convergence des algorithmes d ' adaptation des 
coefficients des filtres adaptatifs etait comparativement 
25 plus faible que la vitesse de convergence de ces algorithmes 
dans le cas de 1' annulation d'echo mono-voie, II a ete 
montre que ce ralentissement de la convergence provient du 
fait que les signaux sonores regus par les haut-parleurs, 
designes par les chiffres de reference 16 1 et 16 2 sur la 
30 figure 2, sont mutuellement corrdl6s . 

Le ralentissement de la convergence entraine 
plusieurs inconvenient s. En particulier, dans un systeme de 
teleconference, les locuteurs situes dans la salle distante 
pergoivent d 1 autant plus longtemps 1 1 echo de leur parole 
35 lors de la mise en route du systeme ou apres une modifica- 
tion acoustique (mouvements des auditeurs par exemple) dans 
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la salle ou est situe l'annuleur d'echo. En outre, en 
annulation d'echo multi-voies, toute modification acoustique 
dans la salle distante perturbe la convergence des filtres 
adaptatifs de l'annuleur d'echo, en raison de la correlation 
5 mutuelle, evoquee precedemment, entre les signaux de parole 
regue, ce qui entraine une reapparition ou une augmentation 
du niveau de 1 1 echo . 

En revanche, on a observe en pratique que lorsque 
des composantes non mutuellement correlees sont prfesentes 

10 sur chaque signal microphonique, elles tendent a accelerer 
la convergence des annuleurs d'echo multi-voies. Un objectif 
de la presente invention est d f exploiter cette propriete des 
composantes mutuellement decorrelees pour ameliorer les 
performances des annuleurs d'echo multi-voies, et ameliorer 

15 par consequent la qualite de la communication dans les 
systemes de transmission de signaux sonores qui mettent en 
oeuvre des annuleurs d'echo multi-voies. Pour cela, le 
principe general de la presente invention consiste a a j outer 
aux voies de signal sonore regu des signaux auxiliaires 

20 mutuellement decorreles, rendus inaudibles par exploitation 
de certaines proprifetes de 1' audition huraaine. 

Plus precisement, la presente invention propose un 
procede d' annulation d'echo sur N voies de signal sonore 
ayant chacune un haut-parleur et un microphone associes, N 

25 etant un entier superieur ou egal a 1, suivant lequel, sur 
chacune des N voies : 

(a) on engendre un signal synthetique ayant les 
caracteristiques spectrales d'un bruit blanc, le spectre de 
ce signal s'etendant sur plusieurs bandes de frequence 

30 adjacentes, et ce signal synthetique etant decorrele des 
signaux synthetiques engendres sur les autres voies ; 

(b) pour chaque bande de frequence, on calcule un 
seuil de masquage frequentiel correspondant au signal 
associe au haut-parleur de la voie consideree en utilisant 

35 des proprietes de la perception auditive humaine ; 

(c) dans chaque bande de frequence, on amene le 



BNSDOCID: <WO 9347276A1_I_> 



WO 98/47276 PCT7FR98/00757 

5 

niveau du signal synthetique a la valeur du seuil de 
masquage frequentiel associd, de fagon a obtenir un signal 
auxiliaire ; 

(d) on additionne le signal auxiliaire au signal 
5 associe au haut-parleur de la voie consideree, les seuils de 
masquage frequentiel ayant ete precfedemment calcules de 
fagon a rendre le signal auxiliaire inaudible, et les 
signaux auxiliaires des N voies etant mutuellement decorre- 
les ; 

10 (e) on fournit le signal obtenu precedemment en 

entree a un -f litre adaptatif dont les coefficients forment 
une estimation de la reponse impulsionnelle du chemin de 
couplage acoustique entre le haut-parleur et le microphone 
associes a la voie de signal sonore consideree ; 

15 (f ) on additionne les signaux obtenus en sortie des 

N filtres adaptatifs respectivement associes a chaque voie 
de microphone et on soustrait le signal resultant du signal 
regu par le microphone associe a cette voie ; 

(g) on calcule une erreur d ? estimation a partir de 
20 la difference obtenue a 1' issue de la soustraction prece- 

dente ; 

(h) on corrige les coefficients des filtres adapta- 
tifs associes a la voie de microphone consideree en fonction 
de 1 ' erreur d ' estimation associee. 

25 Dans un mode particulier de realisation du proc6de, 

pour calculer chaque seuil de masquage frequentiel, 

(bl) on forme des blocs contenant chacun un nombre 
predetermine d r echantillons de signal sonore, deux blocs 
successifs se recouvrant mutuellement sur un nombre prede- 
30 termine d 1 echantillons ; 

(b2) on pondere les echantillons de chaque bloc par 
une f enetre d 1 apodisation ; 

(b3) on calcule la transformee de Fourier de chaque 

bloc ; 

35 (b4) on divise le domaine des frequences en plu- 

sieurs bandes critiques adjacentes ayant des frequences 
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inferieures et superieures determinees ; 

puis, pour chaque bloc et dans chaque bande critique : 

(b5) on additionne l'energie des differentes raies 
spectrales appartenant a la bande critique, de facon a 
obtenir la valeur, pour cette bande critique, d'une fonction 
de repartition de l'energie .; 

(b6) on calcule le produit de convolution de la 
fonction de repartition de l'energie et' d'une fonction 
d'etalement basilaire, obtenue a partir d'une table de 
consultation, de facon a obtenir un spectre etale E ± ; 

(b7) on calcule un indice de tonalite a x a partir des 
energies des raies spectrales appartenant aux differentes 
bandes critiques ; 

(b8) on calcule un facteur de correction 0* a partir 
de 1' indice de tonalite calcule precedemment ; 

(b9) on calcule un seuil de masquage frequentiel T 1 
a partir du spectre etale E, et du facteur de correction 0 ir 
defini comme suit : 

10 log 10 T x = 10 log 10 E L - O t - c t 
ou C ± est un parametre supplementaire de correction de seuil 
calcule a partir de plusieurs caracteristiques des signaux 
sonores des N voies. 

Dans un mode particulier de realisation du procede, 
les caracteristiques servant a calculer le parametre 
supplementaire de correction de seuil C A comprennent le 
niveau respectif des signaux des N voies et un ensemble de 
valeurs de la fonction d ' intercorrelation d'au moins 
certaines paires de ces signaux, le parametre supplementaire 
de correction C ± etant une fonction croissante des valeurs 
de cette fonction d ' intercorrelation. 

En variante, les caracteristiques servant a calculer 
le parametre supplementaire de correction de seuil c x 
peuvent comprendre la somme des energies contenues dans les 
diverses bandes critiques de chaque voie et le rapport des 
energies par bande critique d'au moins certaines paires des 
differentes voies, le parametre supplementaire de correction 
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C± etant d'autant plus important que la somme des energies 
contenues dans les diverses bandes critiques est importante 
et que le rapport des energies par bande critique d'une 
paire de voies est peu different entre les differentes 
5 bandes critiques. 

La presehte invention propose egalement un annuleur 
d 1 echo a N voies de signal sonore ayant chacune un haut- 
parleur et un microphone associes, N etant un entier 
super ieur ou egal a 1, comprenant : 

10 - des moyens de calcul de signaux auxiliaires, 

comportant : 

des moyens de formation de blocs contenant 
chacun un nombre predetermine d T echantillons du 
signal sonore, deux blocs success! fs se recouvrant 
15 mutuellement sur un nombre predetermine d f echantil- 

lons ; 

des moyens de ponderation des echantillons de 
chaque bloc par une fenetre d f apodisation, disposes 
en sortie des moyens de formation de blocs ; 
20 des moyens de calcul de trans formee de 

Fourier, disposes en sortie des moyens de pondera- 
tion ; 

des moyens de calcul des valeurs de 
frequences inferieures et sup6rieures de bandes 
25 critiques adjacentes, ou une table de consultation 

contenant ces valeurs ; 

des moyens pour calculer une fonction de 
repartition de l'energie sur les bandes critiques, 
la valeur, dans une bande critique donnee, de la 
30 fonction de repartition de l'energie etant calcul§e 

en additionnant l'energie des differentes raies 
spectrales appartenant a cette bande critique ; 

une seconde table de consultation, contenant 
des valeurs representant des fonctions d'etalement 
35 basilaire associees chacune a une bande critique 

donnee ; 
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des moyens pour calculer un spectre etal6 E ± , 
en qalculant, pour une bande critique donnee, le 
produit de convolution de la fonction de repartition 
de 1 1 energie et de la fonction d ' etalement associees 
5 a cette bande critique ; 

des moyens de calcul d'un facteur de correc- 
tion 0 A associe a une bande critique donnee ; 

des moyens de calcul d'un seuil de masquage 
frequentiel T L associe a une bande critique donnee, 
10 a partir du spectre etale et du facteur de correc- 

tion Oi associes a cette bande critique, le seuil de 
masquage frequentiel T £ etant defini comme suit : 

10 log 10 T ± = 10 log 10 E ± - 0 ± - C L 
ou C t est un parametre supplementaire de correction 
15 de seuil ; 

des moyens de generation de signaux syntheti- 
ques ayant les caracteristiques spectrales d'un 
bruit blanc ; 

des moyens de controle, connectfes aux N voies 
20 de signal sonore, pour calculer pour chaque bande 

critique le parametre supplementaire de correction 
de seuil C ± a partir de plusieurs caracteristiques 
des signaux sonores des N voies ; 

des moyens pour amener le niveau d'un signal 
25 synthfetique a la valeur du seuil de masquage fre- 

quentiel associe, de fagon a obtenir un signal 
auxiliaire ; 

des moyens de calcul de transformee de Fourier 
inverse ; 

30 des moyens de memoire pour stocker le resultat 

du traitement des blocs successifs ; 

des moyens pour combiner le resultat du 

traitement d 1 un bloc avec le resultat du traitement 

du bloc precedent ; 
35 des moyens de lecture sequent ielle relies 

auxdits moyens de memoire, pour faire passer la 
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cadence des donnees de la cadence bloc a la cadence 
echantillon ; 

N premiers moyens additionneurs respectivement 
disposes sur les N voies de signal sonore, pour 
5 addi*tionner sur chaque voie le signal sonore regu et 

le signal auxiliaire associe ; 

- NxN moyens de filtrage adaptatif , comprenant N 
moyens de filtrage adaptatif associes a chacune des N voies 
de microphone, et ayant chacun des coefficients qui forment 

10 une estimation de la reponse impulsionnelle du chemin de 
couplage acoustique entre le haut-parleur et le microphone 
associes a une des N voies de signal sonore ; 

- des seconds moyens additionneurs, connectes en 
sortie des N moyens de filtrage adaptatif associes a chaque 

15 voie de microphone, pour additionner les signaux de sortie 
de ces N moyens de filtrage adaptatif ; 

- N moyens soustracteurs, disposes chacun en sortie 
des seconds moyens additionneurs et connectes respectivement 
aux microphones des N voies de signal sonore, pour sous- 

20 traire sur chaque voie le signal obtenu en sortie des 
seconds moyens additionneurs du signal regu par le micro- 
phone de cette voie ; 

- N moyens pour calculer sur chaque voie une erreur 
d f estimation a partir du resultat fourni par le moyen 

25 soustracteur connecte au microphone de cette voie ; 

- N moyens pour corriger de fagon iterative les 
coefficients respectifs des N moyens de filtrage adaptatif 
associes a chacune des N voies de microphone, en fonction de 
1' erreur d' estimation associee a chaque voie de microphone. 

30 Les realisations different selon les applications, 

et en particulier, le nombre de haut-parleurs du systeme de 
transmission de signaux sonores considere peut etre diffe- 
rent du nombre de microphones. 

Dans un mode particulier de realisation du disposi- 

35 tif , les caractferistiques servant a calculer le parametre 
supplementaire de correction de seuil C x comprennent le 
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niveau respectif des signaux des N voies et un ensemble de 
valeurs de la fonction d ! intercorrelation d'au moins 
certaines paires de ces signaux, le parametre supplementaire 
de correction C A etant une fonction croissante des valeurs 
5 de cette fonction d ' intercorrelation. 

En variante, les caracteristiques servant a calculer 
le parametre supplementaire de correction de seuil c ± 
peuvent comprendre la somme des energies contenues dans les 
diverses bandes critiques de chaque voie et le rapport des 
10 energies par bande critique d'au moins certaines paires des 

differentes voies, le parametre supplementaire de correction 
Ci etant d'autant plus important que la somme des energies 
contenues dans les diverses bandes critiques est importante 
et que le rapport des energies par bande critique d'une 
paire de voies est peu different entre les differentes 
bandes critiques . 

Dans un mode particulier de realisation du disposi- 
tif, l'annuleur d 1 echo comprend en outre des moyens de 
. calcul d'un indice de tonalite a ± a partir des energies des 
raies spectrales appartenant aux differentes bandes criti- 
ques, lesdits moyens de calcul du facteur* de correction O a 
calculant le facteur de correction 0 ± a partir de 1' indice 
de tonalite a x . 

Dans un mode particulier de realisation du disposi- 
25 tif, les moyens de calcul de 1 1 indice de tonalite fournis- 
sent un indice de tonalite constant, identique pour toutes 
les bandes critiques et def ini par : 

a L = a « minfSFM/SFM^, 1) 
ou min(a, b) designe la plus petite des valeurs a et b, 
ou SFM roax est un parametre de valeur predetermines en dB 
associe a un signal sinusoidal pur, et 

ou SFM = 10 log 10 G/A, ou log 10 designe le logarithme decimal, 
G designe la moyenne geometrique de l'energie sur un nombre 
predetermine de points de la transformee de Fourier et A 
35 designe sa moyenne arithmetique sur le meme nombre de 
points . 
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Dans un mode particulier de realisation, SFM max = 

- 60 dB. 

Dans un mode particulier de realisation, les moyens 
de calcul du facteur de correction 0 L fournissent un facteur 
5 de correction O l defini par : 

0 L = max (SO, a*. (k x + B A ) + (1 - a x ).k 2 ) 
ou max(a, b) designe la plus grande des valeurs a et b, 
ou SO, k x et k 2 sont des parametres de valeur predeterminee 
en dB, 

10 ou a t est 1 ? indice de tonality associe a la bande critique 
consideree, et 

ou B L designe la frequence de la bande critique en Bark. 

Dans un mode particulier de realisation du disposi- 
tif, SO - 24,5 dB, k x = 14,5 dB et k 2 = 5,5 dB. 

15 D'autres caracteristiques et avantages de la 

presente invention apparaitront a la lecture de la descrip- 
tion detaillee qui suit de modes particuliers de realisa- 
tion, donnes a titre d'exemples non limitatifs. La descrip- 
tion se r6fere aux dessins annexes, dans lesquels : 

20 - la figure 1, relative a l'art anterieur et decrite 

en introduction, est une representation sch§matique d'un 
annuleur d'6cho acoustique mono-voie classique ; 

- la figure 2, egalement relative a 1 ' art anterieur 
et decrite en introduction, est une representation schemati- 

25 que partielle d'un annuleur d'echo acoustique stereophonique 
classique, ou une seule voie de microphone est representee ; 

- la figure 3 est un organigramme deer i van t la 
succession generale des etapes du procede de la presente 
invention dans un mode particulier de realisation ; 

30 - la figure 4 represente de fagon schematique le 

dispositif de la presente invention dans un mode particulier 
de realisation ; 

- la figure 5 represente de fagon schdmatique la 
structure du module de calcul de signaux auxiliaires compris 

35 dans le dispositif de la presente invention, dans un mode 
particulier de realisation ; 
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- la figure 6a est un graphe illustrant un premier 
exemple d' evolution au cours du temps de 1 ' ecart relatif 
entre les reponses impulsionnelles reelles des chemins de 
couplage acoustique et les reponses impulsionnelles estimees 

5 par les filtres adaptatifs d'un annuleur d'^cho stereophoni- 

que conforme a un mode particulier de realisation de la 
presente invention ; 

- la figure 6b est un graphe analogue a celui de la 
figure 6a , etabli pour un annuleur d'echo stereophonique 

10 classique ; 

- la figure 7a est un graphe illustrant un deuxieme 
exemple d' evolution au cours du temps de 1 1 ecart relatif 
entre les reponses impulsionnelles reelles des chemins de 
couplage acoustique et les reponses impulsionnelles estimees 

15 par les filtres adaptatifs d'un annuleur d'echo stereophoni- 

que conforme a un mode particulier de realisation de la 
presente invention ; 

- la figure 7b est un graphe analogue a celui de la 
figure 7a, etabli pour un annuleur d'echo stereophonique 

20 classique. 

Le procede et le dispositif de la presente invention 
s'appliquent a un signal sonore quelconque. A titre d 1 exem- 
ple non limitatif, il peut s'agir aussi bien d'un signal de 
parole que de musique. 

25 0n va d'abord decrire 1 1 organisation generale du 

procede d' ambulation d'echo multi-voies dans un mode 
particulier de realisation, a l'aide de la figure 3. 

On considere N voies de signal sonore dont chacune 
comporte un haut-parleur et un microphone. 

30 Sur chaque voie de signal sonore regu, on synthetise 

tout d'abord un signal ayant les caracteristiques spectrales 
d'un bruit blanc et dont le spectre s'etend sur plusieurs 
. bandes de frequence adjacentes, les signaux synthetiques des 
N voies etant mutuellement decorreles. 
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Ensuite, pour chaque bande de frequence, on calcule 
un seuil de masquage frequentiel correspondant au signal 
sonore regu par le haut-parleur de la voie consideree. 

On rappelle que l'effet de masquage consiste en ce 
5 qu'un son en presence d'un autre peut devenir partiellement 
ou completement inaudible. Lorsque ce phenomene se produit 
dans le domaine frequentiel, on parle de masquage simultane 
ou masquage frequentiel. II existe un seuil, appele seuil de 
masquage frequentiel, en dega duquel toutes les composantes 

10 frequentielles d'un son masque sont inaudibles. Pour plus de 
detail sur le phenomene de masquage, on se reportera 
utilement a l'ouvrage de E. ZWICKER et R. FELDTKELLER 
intitule "Das Ohr als Nachrichtenempfanger", Stuttgart, West 
Germany, Hirzel Verlag, 1967. 

15 Dans un mode particulier de realisation, le calcul 

des seuils de masquage frequentiel s'effectue comme suit. 
Tout d'abord, on transforme la cadence du flot de donn6es de 
signal sonore en passant d'une cadence d 1 echanti lions , par 
exemple de 16 kHz, a une cadence de blocs, par exemple 

20 chacun d'une duree de 32 ms. Chaque bloc contient done un 
nombre predetermine d 1 echantillons . Deux blocs successifs se 
recouvrent mutuellement sur un nombre predetermine d T echan- 
tillons, par exemple sur une duree de 16 ms . 

Ensuite, pour ameliorer l f estimation ulterieure du 

25 spectre, on pond<lre les echantillons de chaque bloc par une 
fenetre d 1 apodisation, par exemple du type fenetre de 
Hamming ou de Harming. 

Puis on calcule la transformee de Fourier de chaque 
bloc sur un- nombre predetermine de points, par exemple 512 

30 points pour des blocs de 32 ms. On peut utiliser a cet effet 
un algorithme de transformee de Fourier rapide. 

On divise ensuite le domaine des frequences auquel 
le signal sonore regu est susceptible d 1 appartenir en 
plusieurs bandes de frequences adjacentes ayant des frequen- 

35 ces inferieures et superieures determinees. Ces bandes de 
frequences sont appelees bandes critiques. Le tableau ci- 
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dessous fournit, a titre d'exemple non limitatif, la 
decomposition du domaine des frequences de la parole en 
bandes critiques (bornes inferieures et superieures et 
frequences centrales) etablie par E. ZWICKER et R. FELDTKEL- 
LER et figurant dans l'ouvrage precite. La largeur de chaque 
bande critique est par convention egale a 1 Bark. Le tableau 
donne done la correspondance entre l'echelle frequentielle 
des Hertz et celle des Bark. 



15 
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25 



30 



35 



40 



Bande 
critique 



1 
2 
3 
4 
5 
6 
7 
8 
9 

10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 



Borne 
inferieure 
(Hz) 



0 
100 
200 
300 
400 
510 
630 
770 
920 
1080 
1270 
1480 
1720 
2000 
2320 
2700 
3150 
3700 
4400 
5300 
6400 
7700 
9500 
12000 



Centre 
(Hz) 



50 
150 
250 
350 
450 
570 
700 
840 
1000 
1170 
1370 
1600 
1850 
2150 
2500 
2900 
3400 
4000 
4800 
5800 
7000 
8500 
10500 
13500 



Borne 
superieure 
(Hz) 



100 
200 
300 
400 
510 
630 
770 
920 
1080 
1270 
1480 
1720 
2000 
2320 
2700 
3150 
3700 
4400 
5300 
6400 
7700 
9500 
12000 
15500 



Au lieu d'appliquer une transformee de Fourier comme 
decrit ci-dessus, on peut, en variante, utiliser 
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decomposition par bancs de filtres, ce qui permet d ' obtenir 
directement la decomposition en bandes critiques. 

L'etape suivante consiste a additionner l'energie 
des differentes raies spectrales appartenant a chaque bande 
5 critique, de fagon a obtenir la valeur, dans cette bande 
critique, d'une fonction de. repartition de l'energie. 

On calcule ensuite le spectre etale E x , d6fini comme 
le produit de convolution de la fonction de repartition de 
1 1 energie par bande critique avec la fonction dite d 1 etale- 

10 ment basilaire. La fonction d'etalement basilaire E b 
correspond a la reponse f requentielle en energie de la 
membrane basilaire de l f oreille humaine. Elle represente 
1 ! excitation resultant de la perception d 1 une tonale, c'est- 
a-dire d'un son de frequence pure, et peut etre definie par 

15 l f equation * suivante, figurant egalement dans l'ouvrage 
precite de E. ZWICKER et R. FELDTKELLER : 
10.1og 10 E b = 15;81+7,5.(B+0 / 474)-17,5.(l+(B+0,474) 2 ) 1/2 
ou log 10 represente le logarithme decimal et B est la 
frequence en Bark. Dans le mode particulier de realisation 

20 decrit, les valeurs de cette fonction sont lues dans une 
table de consultation. 

On a montre exp6rimentalement qu 1 il est plus 
difficile de masquer un bruit par une tonale que 1' inverse. 
Pour tenir compte de cette asymfetrie des seuils de masquage, 

25 on applique a chaque seuil un facteur de correction 0 i 
associe a une bande critique donnee B 1 et qui depend de 
l f aspect tonal ou non tonal du signal sonore. Pour evaluer 
cet aspect tonal ou non tonal, on utilise, dans un mode 
particulier de realisation, une mesure notee SFM (en anglais 

30 Spectral Flatness Measure), definie comme le rapport entre 
la moyenne geometrique de l'energie sur un nombre predeter- 
mine de points de la transformer de Fourier et sa moyenne 
arithmetique sur le meme nombre de points, par exemple 256 
points pour une trans formee de Fourier calculee sur 512 

35 points. 
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La valeur SFM obtenue est comparee a celle d ' un 
signal sinusoidal pur, pour lequel, a titre d'exemple non 
limitatif , on peut poser par convention SFM = SFM^ = 
- 60 dB. 

5 On calcule alors un indice de tonalite a if qui, dans 

un mode particulier de realisation, est constant et identi- 
que pour toutes les bandes critiques et defini par : 

a A - a = min ( SFM/SFM^, 1) 
ou min (a, b) designe la plus petite des valeurs a et b. 
10 L' indice a ± est compris entre 0 et 1 et est plus eleve pour 
les composantes tonales du signal que pour les composantes 
non tonales. Toute autre methode de calcul de 1 ' indice de 
tonalite jugee appropriee est envisageable. 

Le facteur de correction O t associe a la bande 
15 critique de frequence B 4 en Bark est defini a partir de 
1 * indice de tonalite a L comme suit : 

Oi = max (SO, ^.(^ + B 4 ) + (1 - a t ) .k 2 ) 
ou max (a, b) designe la plus grande des valeurs a et b et 
ou SO, k x et k 2 sont des parametres de valeur predeterminee 
en dB. Typiquement, on peut choisir SO = 24,5 dB, kj = 14,5 
dB et k 2 = 5,5 dB . 

Le facteur de correction Oi est une fonction crois- 
sante du rang de la bande critique dans le cas de composan- 
tes tonales dominantes et accorde plus de poids a celles-ci 
25 qu ' aux composantes non tonales . 

On applique en outre aux seuils de masquage un 
parametre supplementaire de correction C x calcule a partir 
de plusieurs caracteristiques des signaux sonores des N 
voies. Le parametre C t peut etre commun a tous les seuils de 
masquage correspondant aux differentes bandes critiques, ou 
bien avoir une valeur specif ique pour chaque bande critique. 

Dans un mode particulier de realisation, on calcule 
le parametre supplementaire de correction de seuil C 4 a 
partir des niveaux des signaux sonores recus sur les N voies 
35 et d'un ensemble de valeurs de la fonction d ' intercorrela- 
tion d'au moins certaines paires de ces signaux. Dans ce 
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mode particulier de realisation, le parametre supplementaire 
de correction de seuil C L est une f onction croissante des 
valeurs de la f onction d 1 intercorrelation. 

A titre d'exemple non limitatif, dans une installa- 
5 tion pratique typique de teleconference a plusieurs micro- 
phones et plusieurs haut-parleurs, les microphones sont 
disposes de telle sorte que la correlation entre signaux 
captes est importante entre deux microphones adjacents et 
faible entre microphones distants. Cette correlation 

10 variable se retrouve, moyennant des modifications eventuel- 
les dues aux traitements durant la transmission, a la 
reception, entre les voies adjacentes d'une part, a savoir, 
entre deux haut-parleurs directement voisins, et entre les 
voies distantes d' autre part, a savoir, entre les autres 

15 paires de haut-parleurs distants. 

Ainsi, on peut se limiter a utiliser, pour le calcul 
des parametres C if les caracteristiques de niveau et 
d 1 intercorrelation des paires de voies correspondant aux 
paires de haut-parleurs directement voisins et aux deux 

20 paires adjacentes (si ces dernieres existent, compte tenu du 
nombre total de haut-parleurs). Ainsi, si on considere une 
paire de haut-parleurs directement voisins (HP k , HP k+1 ), on 
utilisera leurs caracteristiques de niveau et d 1 intercorre- 
lation ainsi que celles des paires adjacentes (HP k-1 , HP k ) 

25 (HP k „! et HP k 6tant directement voisins) et (HP k4l , HP k . 2 ) (HP k ^ 

et HP k+2 etant directement voisins). 

Dans un autre mode de realisation, on calcule le 
parametre supplementaire de correction Ci a partir de la 
somme des energies contenues dans les diverses bandes 

30 critiques de chaque voie et le rapport des energies par 
bande critique d'au moins certaines paires des differentes 
voies . 

De fagon analogue au mode de realisation precedent, 
on peut se limiter a utiliser les caracteristiques de somme 
35 et de rapport d'fenergie des paires de voies correspondant 
aux paires de haut-parleurs directement voisins et aux deux 
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paires adjacentes. Le parametre C ± est d'autant plus 
impor-tant que la somme des energies contenues dans les 
diverses bandes critiques est importante et que le rapport 
des energies par bande critique d 1 une paire de voies est peu 
5 different entre les differentes bandes critiques. 

Ainsi, lorsque 1 1 intercorrelation entre les signaux 
sonores de deux voies donnees est forte, et/ou que ces deux 
signaux ont des energies comparables, le parametre C L est 
eleve, c'est-a-dire qu'une reduction importante du seuil est 

10 appliquee, typiquement de 6 a 10 dB. A 1 ? inverse, lorsque 
les signaux sont peu correles et/ou sont d f energies sensi- 
blement differentes, une reduction moins importante du seuil 
est appliquee, typiquement de 0 a 6 dB. Cette correction 
tient compte des proprietes de 1' audition. Elle a pour but 

15 de rendre inaudibles les signaux auxiliaires mentionnes en 

introduction et dont I'obtention est detaillee ci-apres, 
tout en permettant de leur donner le niveau le plus eleve 
possible. 

Pour chaque bloc et dans chaque bande critique, on 
obtient le seuil de masquage frequentiel T x a partir du 
spectre etale E ± et du facteur de correction 0 ± definis 
precedemment , et du parametre supplementaire de correction 
de seuil C if comme suit : 

10 log 10 T ± . 10 log 10 E t - 0 ± - c ± 
25 Comme le montre la figure 3, l'etape suivante du 

procede consiste a obtenir, sur chaque voie, un signal 
auxiliaire. Pour cela, dans un mode particulier de realisa- 
tion, on amene dans chaque bande critique le niveau du 
signal synthetique a la valeur du seuil de masquage frequen- 
30 tiel associe. 

On calcule ensuite la trans formee de Fourier inverse 
correspondant au resultat du traitement du bloc courant, par 
exemple au moyen d'un algorithme de transformee de Fourier 
inverse rapide. Dans le cas ou on n»a pas utilise une 
35 transformee de Fourier, mais une decomposition par bancs de 
filtres, comme decrit precedemment, on effectue, au lieu 
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d 1 une -transformation de Fourier inverse, une etape de 
synthese correspondant a une transformation inverse de la 
decomposition par bancs de filtres. 

On combine ensuite le resultat du traitement des 
5 blocs precedents avec le resultat du traitement du bloc 

courant, afin de reconstruire le bloc de bruit masque 
courant, par exemple au moyen d'un procede connu de synthese 
par recouvrement et addition (en anglais " over 1 ap- add " ) ou 
par recouvrement et sauvegarde (en anglais " overlap- save" ) . 
10 On repasse ensuite de la cadence de blocs a la 

cadence d f echantillons . 

L f ensemble des traitements decrits precfedemment peut 
etre effectue par un processeur de traitement de signal en 
temps reel. 

15 Dans un mode particulier de realisation, 1' etape 

suivante du .procede consiste, sur chaque voie, a additionner 
le signal auxiliaire obtenu au signal associe au haut- 
parleur de cette voie, comme 1 ■ indique la figure 3. Par 
construction, les signaux auxiliaires des N voies sont 

20 mutuellement decorreles et, du fait que les signaux auxi- 
liaires ne depassent pas les seuils de masquage frequentiel 
associes, ils sont inaudibles. 

Dans le mode de realisation decrit, chaque signal 
auxiliaire est additionne au signal de haut-parleur associe 

25 apres avoir ete reconstitue dans sa totalite a partir des 
differentes bandes critiques. En variante, cette addition 
peut tout aussi bien etre effectuee dans chacune des 
differentes bandes critiques. 

On effectue ensuite une etape globale de filtrage 

30 adaptatif, comme 1 1 indique la figure 3. Cette etape est 
detaillee ci-dessous dans un mode particulier de realisa- 
tion. 

On fournit tout d'abord la somme du signal auxi- 
liaire et du signal de haut-parleur en entree a un filtre 
35 adaptatif dont les coefficients forment une estimation de la 
reponse impulsionnelle du chemin de couplage acoustique 
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entre le haut-parleur et le microphone associes a la voie de 
signal sonore consideree. 

On procede ainsi sur chaque voie, puis on additionne 
les signaux obtenus en sortie des N filtres adaptatifs 
5 associes a chaque voie de microphone. On soustrait ensuite 

le signal resultant, sur chaque voie de microphone, du 
signal regu par le microphone associe a cette voie. 

A partir du signal de difference ainsi obtenu, on 
calcule sur chaque voie de microphone une erreur d'estima- 
10 tion et on corrige les coefficients des filtres adaptatifs 

associes a cette voie de microphone en fonction de cette 
erreur d 1 estimation. 

L'etape de filtrage adaptatif detaillee ci-dessus 
est realisee dans le domaine temporel. En variante, on peut 
15 realiser le filtrage adaptatif par blocs et non par echan- 

tillons, en particulier en effectuant 1' ensemble des traite- 
ments correspondents dans le domaine de Fourier. Pour une 
description de procedes de filtrage adaptatif frequentiel, 
on se reportera utilement a 1' article de J. PRADO et E. 
20 MOULINES intitule "Frequency domain adaptative filtering 

with applications to acoustic echo cancellation", Annales 
des Telecommunications, 49, n° 7-8, 1994, pages 414-428. 

On decrit dans ce qui suit un mode particulier de 
realisation d'un annuleur d'echo conforme a la presente 
25 invention, a l'aide des figures 4 et 5. 

Dans ce mode particulier de realisation, le nombre 
de haut-parleurs est identique au nombre de microphones. 
Cependant, diverses realisations sont possibles selon les 
applications, et en particulier, le nombre de haut-parleurs 
30 peut etre different du nombre de microphones. 

Comme le montre la figure 4, 1' annuleur d'echo est 
associe a N haut-parleurs 30 a , . . ., 30 N et N microphones 32 x , 
32 N . II comprend un module 34 de calcul de signaux 
auxiliaires dont la structure detaillee sera decrite a 
35 l f aide de la figure 5. De fagon globale, le module 34 de 

calcul de signaux auxiliaires comprend un module de controle 
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36 connecte aux N voies de haut-parleur, N modules 38!, 

38 N de production de signaux auxiliaires respectivement 

connectes aux N voies de haut-parleur et au module de 

controle 36, et N premiers modules additionneurs 40^ , 

40 N respectivement disposes sur les N voies de haut-parleur 
et recevant chacun en entree le signal sonore regu sur la 
voie consideree et le signal de sortie du module de produc- 
tion de signaux auxiliaires associe a cette voie. 

L'annuleur d'echo represents sur la figure 4 
comprend en outre N x N modules de filtrage adaptatif 42 A 4 
ou i et j sont des entiers variant de 1 a N, les N modules 

de filtrage adaptatif 42^ 42 i>H etant associes a la 

ieme voie de microphone, 1 s i < N, et etant respectivement 
associes aux N voies de haut-parleur, comme le montre la 
figure 4. Pour la ieme voie de microphone, 1 s i < n, les 
sorties des N modules de filtrage adaptatif 42 x lf ... 42 
sont connectees a un second module additionneur 44 x . La 
sortie de chaque second module additionneur 44 x , 1 < i s n, 
est connectee a 1' entree soustractrice d'un module soustrac- 
teur 46 A recevant sur son entree non soustractrice le signal 
recu par le ieme microphone 32*. La sortie de chaque module 
soustracteur 46 x commande 1' ensemble des N modules de 

filtrage adaptatif 42 1>x 42 1>N , par 1 • intermediate d'un 

module (non represents ) permettant de calculer une erreur 
d' estimation a partir du resultat fourni par le module 
soustracteur 46^ et par 1 ' intermediaire d'un module (non 
represents) pour corriger de facon iterative les coeffi- 
cients des N modules de filtrage adaptatif 42 t x , , 42 ± „, 

en fonction de 1' erreur d' estimation obtenue. 

La figure 5 represente plus en detail la structure 
du sous-ensemble du module 34 de calcul de signaux auxiliai- 
res forme par le module de contr61e 36 et les N modules 38 x , 
38 N de production de signaux auxiliaires. 

Le signal sonore regu sur une voie de haut-parleur 
donnee est fourni en entree a un module 48 de formation de 
blocs. Le module 48 fournit en sortie des blocs contenant 
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chacun un nombre predetermine d'echantillons du signal 
sonore regu.^ Deux blocs successifs se recouvrent mutuelle- 
ment sur un nombre predetermine d'echantillons. Les exemples 
numeriques donnes dans la description d 1 un mode particulier 
5 de realisation du procede restent valables. 

Le module 48 de formation de blocs est connecte a un 
module de ponderation 50 qui applique a chaque bloc une 
fenetre d ' apodisation, par exemple du type fenetre de 
Hamming ou de Hanning. 

10 Le module de ponderation 50 est connecte a un module 

52 de calcul de transformee de Fourier. A titre d ' exemple 
non limitatif , le module 52 peut effectuer une transformee 
de Fourier rapide sur 512 points. 

Le module 34 de calcul de signaux auxiliaires 

15 comprend egalement un module 54 de calcul des valeurs des 

frequences inferieures et superieures de bandes critiques 
adjacentes couvrant le domaine de frequences auquel appar- 
tient le signal sonore regu. En variante, le module 54 peut 
consister en une table de consultation contenant ces 

20 valeurs. 

Le module 54 est connecte a un module 56 de calcul 
de fonctions de repartition de l'energie. L ' entree du module 
56 est egalement connectee au module 52 de calcul de 
transformee de Fourier. Le module 56 calcule et additionne 

25 l'energie des differentes raies spectrales appartenant a 
chaque bande critique en utilisant les valeurs des frequen- 
ces limites des bandes critiques fournies par le module 54. 
Le module 56 fournit en sortie, pour chaque bande critique, 
la valeur d'une fonction de repartition de l'energie. 

30 Une seconde table de consultation 58 contient des 

valeurs representant des fonctions d'etalement basilaire 
associees chacune a une bande critique donnee. La definition 
de ces fonctions a ete donnee plus haut en relation avec un 
mode particulier de realisation du procede. 

35 La seconde table de consultation 58 est connectee a 

I 1 entree d'un module 60 permettant de calculer un spectre 
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etal6 Ei. L 1 entree du module 60 est egalement connectee au 
module 56 de calcul de f ©notions de repartition de l'ener- 
gie. Le module 60 calcule, pour une bande critique donnee, 
le produit de convolution de la fonction de repartition de 
5 l'energie et de la fonction d'etalement, de fagon a fournir 
en sortie le spectre etale pour cette bande critique. 

Dans le mode particulier de realisation represents 
sur la figure 5, le module 34 de calcul de signaux auxiliai- 
res comprend en outre un module 62 de calcul d T indices de 

10 tonalite a t respectivement associes a des bandes critiques 
donnees. Le module 62 est connecte en sortie du module 52 de 
calcul de transformee de Fourier et utilise les valeurs des 
energies des raies spectrales appartenant aux differentes 
bandes critiques . 

15 Dans ce mode particulier de realisation, le module 

62 de calcul d' indices de tonalite est connecte a un module 
64 de calcul de facteurs de correction O a respectivement 
associes a* des bandes critiques donnees. Le module 64 
calcule les facteurs de correction O t a partir des indices 

20 de tonalite a ± . 

Dans un mode particulier de realisation, le module 
62 de calcul d 1 indices de tonalite fournit un indice de 
tonalite a t constant, identique pour toutes les bandes 
critiques, defini par : 

25 a L . - a - min( SFM/SFM max , 1) 

ou min(a, b) designe la plus petite des valeurs a et b, 

ou SFM„ ax est un parametre de valeur predeterminee en dB 

associe a un signal sinusoidal pur, et 

ou SFM = 10 log 10 G/A, ou log 10 designe le logarithme decimal, 
30 G designe la moyenne geom^trique de l'energie sur un nombre 
predetermine de points de la transformee de Fourier et A 
designe sa moyenne arithmetique sur le meme nombre de 
points. Dans ce mode de realisation, le module 64 de calcul 
de facteurs de correction O t fournit en sortie pour chaque 
35 bande critique un facteur de correction 0 i dfefini par : 

0 L = max (SO, a L . (k x + B t ) + (1 - a A ).k 2 ) 
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ou max(a, b) designe la plus grande des valeurs a et b, 

ou SO, k a et. k 2 sont des parametres de valeur predeterminee 

en dB, 

ou a x est: l'indice de tonalite associe a la bande critique 
5 consideree, et 

ou B ± designe la frequence de la bande critique en Bark. 
A titre d'exemple non limitatif, les parametres SFM^, SO, 
k x et k 2 peuvent avoir les valeurs numeriques indiguees plus 
haut en relation avec un mode particulier de realisation du 
10 procede. 

Le module 64 de calcul de facteurs de correction O t 
et le module 60 de calcul de spectres etales sont connectes 
a I 1 entree d'un module 66 de calcul de seuils de masquage 
frequentiel T ± respectivement associes a des bandes criti- 
15 ques donnees. Le module 66 fournit en sortie pour chaque 

bande critique un seuil de masquage frequentiel defini 
comme suit : 

10 log 10 Ti - 10 log 10 E t - 0 ± - 
ou Ei est le spectre etale et 0 L le facteur de correction 

20 associes a cette bande critique, et ou Ci est un parametre 

supplementaire de correction de seuil. 

Le parametre C ± est fourni par le module de controle 
36 mentionne plus haut. Les connexions du module 36 aux N 
voies de haut-parleur, representees sur la figure 4, ont et6 

25 omises sur la figure 5 pour plus de clarte. 

Dans un mode particulier de realisation, le module 
de controle 36 calcule le parametre supplementaire de 
correction de seuil C A pour chaque bande critique a partir 
du niveau respectif des signaux des N voies et d f un ensemble 

30 de valeurs de la fonction d 1 intercorrelation d'au moins 
certaines paires de ces signaux, comme on l'a decrit plus 
haut en relation avec un mode particulier de realisation du 
procede. On a vu que dans ce cas, le parametre C L est une 
fonction croissante des valeurs de la fonction d ■ intercorre- 

35 lation. 
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Dans un autre mode particulier de realisation, le 
module de control e 36 calcule le parametre C A pour chaque 
bande critique a partir de la somme des energies contenues 
dans les diverses bandes critiques de chaque voie et du 
rapport des energies par bande critique d 1 au moins certaines 
paires des differentes voies, comme on 1 T a egalement decrit 
en relation avec un mode particulier de realisation du 
procede. On a vu que dans ce cas, le parametre C L est 
d'autant plus important que la somme des energies contenues 
dans les diverses bandes critiques est importante et que le 
rapport des energies par bande critique d'une paire de voies 
est peu different entre les differentes bandes critiques. 

Le module 34 de calcul de signaux auxiliaires 
comporte de plus un module 68 de generation de signaux 
synthetiques . Le module produit, pour chaque voie de signal 
sonore, un signal ayant les caracteristiques spectrales d'un 
bruit blanc et decorrel6 des signaux synthetiques produits 
pour les autres voies. 

La sortie du module 66 de calcul de seuils de 
masquage f requentiel est connectee a un module de conversion 
frequentielle 69. Le module 69 effectue la conversion des 
seuils de masquage fournis par le module 66 de Bark en 
Hertz . 

Un module 70 regoit en entree les signaux syntheti- 
ques produits par le module 68 et les seuils de masquage 
frequentiel calcules par le module 66 et convertis par le 
module 69. Le module 70 amene le niveau de chaque signal 
synthetique a la valeur du seuil de masquage frequentiel 
associe, de fagon a obtenir un signal auxiliaire. 

La sortie du module 70 est connectee a l 1 entree d'un 
module 72 de calcul de transformee de Fourier inverse. 

La sortie du module 72 est connectee a l 1 entree d'un 
module de synthese 74. Le module de synthese 74 est connecte 
& une m6moire vive (non representee) qui stocke le resultat 
du traitement des blocs d ' echantillons successifs. Le module 
de synthese 74 combine le resultat du traitement d'un bloc 
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avec le resultat du traitement du bloc precedent: fourni par 
la memoire vive, par exemple en mettant en oeuvre un procede 
de synthese par recouvrement et addition (en anglais 
"overlap-add") ou recouvrement et sauvegarde (en anglais 
5 "overlap-save" ) . 

Un module 76 de lecture sequentielle est connecte au 
module de synthese 74 et a la memoire vive, et assure le 
passage de la cadence des donnees de la cadence de blocs a 
une cadence d ' echantillons . 

10 L ' ensemble des modules 48 a 76 qui viennent d'etre 

decrits est equivalent a la serie de modules 38 lr . . . , 38 N de 
production de signaux auxiliaires representes de fagon 
simplifiee sur la figure 4. En effet, chacun des modules 48 
& 76 contribue a produire un signal auxiliaire pour chaque 

15 voie de haut-parleur. 

Ainsi, la presente invention permet d'injecter a 
chaque instant sur chaque voie de signal sonore regu un 
signal auxiliaire de niveau maximal, assurant efficacement 
1 ' accroissement de la vitesse de convergence, et dont le 

20 niveau est rendu par construction inferieur au seuil de 

masquage par le signal sonore regu. 

Le procede et le dispositif de la presente invention 
s'appliquent egalement en annulation d'echo mono-voie (cas 
N = 1), pour ameliorer la convergence du filtre adaptatif. 

25 Us sont applicables a tous les types d f annulation 

d'echo acoustique, en particulier aussi bien dans le domaine 
temporel que dans le domaine freguentiel, ou encore en sous- 
bandes . 

Les figures 6a, 6b et 7a, 7b illustrent deux 
30 exemples de resultats comparatifs obteniis a l'aide du 
procede et du dispositif d 1 annulation d'echo conformes a un 
mode particulier de realisation de la presente invention. 

Les graphes des figures 6a a 7b comprennent tous 
deux courbes. Les courbes en tirets concernent la voie 
35 ster<§ophonique droite et les courbes en pointilles concer- 
nent la voie stereophonique gauche. L'algorithme utilise 
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pour 1' adaptation des coefficients des filtres adaptatifs 
est l'algorithme du gradient stochastique normalise (NLMS). 
Les filtres adaptatifs ont une taille de 1024 coefficients. 
L f axe des abscisses represente le nombre d 1 iterations de 
5 l'algorithme d 1 adaptation, divise par un facteur 256. L'axe 
des ordonnees represente l'ecart relatif, en dB, entre les 
reponses impulsionnelles reelles des chemins de couplage 
acoustique et les reponses impulsionnelles des filtres 
adaptatifs d'un annuleur d'echo stereophonique. 

10 Le premier exemple comparatif, illustre par les 

figures 6a et 6b, correspond a un locuteur masculin. Le pas 
d' adaptation de l'algorithme NLMS est de 0,25. La figure 6a 
montre l'ecart relatif obtenu en utilisant la presente 
invention et la figure 6b montre l'ecart relatif obtenu en 

15 utilisant un procede et un dispositif classiques d f ambula- 
tion d'echo, ne faisant notamment pas intervenir de signaux 
auxiliaires. 

Le deuxieme exemple comparatif, illustre par les 
figures 7a et 7b, correspond a un locuteur feminin. Le pas 

20 d' adaptation de l'algorithme NLMS est de 0,125. Les figures 
7a et 7b sont les analogues respectifs des figures 6a et 6b. 

Dans les deux exemples comparatif s, 1 1 utilisation de 
la presente invention assure tou jours un ecart relatif 
inferieur a 1 ( c 1 est-a-dire negatif en dB), ce qui n'est pas 

25 le cas du procede et du dispositif classiques pour la voie 
droite (courbes en tirets ) . De plus, dans les cas d* utilisa- 
tion de la presente invention, l'ecart sur la voie gauche 
(courbes en pointilles ) est inferieur de plusieurs dB a 
l'ecart observe sur les courbes correspondant au procede et 

30 au dispositif classiques. La presente invention presente 
done des performances sup6rieures, permettant de garantir un 
bon comportement de poursuite de 1 ' annuleur d'echo en 
presence de variations des chemins de couplage acoustique. 

Par ailleurs, dans les deux exemples, la correction 

35 appliquee aiix seuils de masquage est de 10 dB. On pourrait 
tol6rer une reduction de quelques dB de cette correction, ce 
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qui ameliorerait encore les performances du procede et du 
dispositif proposes, alors que la parole pergue, bien que 
legerement modifiee par le bruit: ajoute, conserverait une 
qualite "tout a fait acceptable dans une application du type 
5 teleconference. 
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RE VEND I CAT I ONS 

1. Proc6de d ' annulation d'echo s.ur N voles de signal 
sonore ay ant chacune un haut-parleur et un microphone 
5 associes, N. etant un entier superieur ou egal & 1, suivant 
lequel, sur chacune des N voies : 

(a) on engendre un signal synthetique ayant les 
caracteristiques spectrales d'un bruit blanc, le spectre de 
ce signal s f etendant sur plusieurs bandes de frequence 

10 adjacentes, et ce signal synthetique etant decorrele des 
signaux synthetiques engendres sur les autres voies ; 

(b) pour chaque bande de frequence, on calcule un 
seuil de masquage f requentiel correspondant au signal 
associe au haut-parleur de la voie consideree en utilisant 

15 des proprietes de la perception auditive humaine ; 

(c) dans chaque bande de frequence, on amene le 
niveau du signal synthetique a la valeur du seuil de 
masquage f requentiel associe, de fagon a obtenir un signal 
auxiliaire ; 

20 (d) on additionne le signal auxiliaire au signal 

associe au haut-parleur de la voie consideree, les seuils de 
masquage f requentiel ayant ete precedemment calcules de 
fagon a rendre le signal auxiliaire inaudible, et les 
signaux auxiliaires des N voies etant mutuellement decorre- 

25 les ; 

(e) on fournit le signal obtenu precedemment en 
entree a un filtre adaptatif dont les coefficients forment 
une estimation de la reponse impulsionnelle du chemin de 
couplage acoustique entre le haut-parleur et le microphone 

30 associes a la voie de signal sonore consideree ; 

(f ) on additionne les signaux obtenus en sortie des 
N f litres adaptatif s associes a chaque voie de microphone et 
on soustrait le signal resultant du signal regu par le 
microphone assocife a cette voie ; 
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(g) on calcule une erreur d 1 estimation a partir de 
la difference obtenue a l f issue de la soustraction prece- 
dente ; 

(h) on corrige les coefficients des filtres adapta- 
5 tifs associes a la voie de microphone consideree en fonction 

de 1 1 erreur d 1 estimation associee . 

2. Procede selon la revendication 1, caracterise en 
ce que, pour calculer chaque seuil de masquage fr6quentiel, 
(bl) on forme des blocs contenant chacun un nombre 
10 predetermine d ' echantillons de signal sonore, deux blocs 
success! fs se recouvrant mutuellement sur un nombre prede- 
termine d 1 echantillons ; 

(b2) on pondere les echantillons de chaque bloc par 
une fenetre d ' apodisation ; 
15 (b3) on calcule la transformee de Fourier de chaque 

bloc ; 

(b4) on divise le domaine des frequences en plu- 
sieurs bandes critiques adjacentes ayant des frequences 
inferieures et superieures determinees ; 
20 puis, pour chaque bloc et dans chaque bande critique : 

(b5) on additionne 1 ' energie des differentes raies 
spectrales appartenant a la bande critique, de fagon a 
obtenir la valeur, pour cette bande critique, d'une fonction 
de repartition de 1 1 energie ; 
25 (b6) on calcule le produit de convolution de la 

fonction de repartition de 1' energie et d'une fonction 
d'etalement basilaire, obtenue a partir d'une table de 
consultation, de fagon a obtenir un spectre etale E A ; 

(b7) on calcule un indice de tonalite a x a partir des 
30 energies des raies spectrales appartenant aux differentes 
bandes critiques ; 

(b8) on calcule un facteur de correction 0 L a partir 
de 1 ' indice de tonalite calcule precedemment ; 

(b9) on calcule un seuil de masquage frequentiel T A 
35 a partir du spectre etale E t et du facteur de correction 0^ 
defini comme suit : 
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10 log 10 T L = 10 log 10 E A - O a - C A 
ou C x est un.parametre supplementaire de correction de seuil 
calcule a partir de plusieurs caracteristiques des signaux 
sonores des N voies. 
5 3. Procede selon la revendication 1 ou 2, caracte- 

ris6 en ce que les caracteristiques servant & calculer ledit 
parametre supplementaire de correction de seuil C L compren- 
nent le niveau respect if des signaux des N voies et un 
ensemble de valeurs de la fonction d 1 intercorrelation d'au 

10 moins certaines paires de ces signaux, ledit parametre 
supplementaire de correction C t etant une fonction crois- 
sante des valeurs de cette fonction d ' intercorrelation. 

4. Procede selon la revendication 1 ou 2, caract6- 
rise en ce que les caracteristiques servant a calculer ledit 

15 parametre supplementaire de correction de seuil C t compren- 
nent la somme des energies contenues dans les diverses 
bandes critiques de chaque voie et le rapport des energies 
par bande critique d'au moins certaines paires des differen- 
tes voies, ledit parametre supplementaire de correction C L 

20 etant d 1 autant plus important que la somme des energies 
contenues dans les diverses bandes critiques est importante 
et que le rapport des energies par bande critique d'une 
paire de voies est peu different entre les differentes 
bandes critiques. 

25 5. Annuleur d'echo a N voies de signal sonore ayant 

chacune un haut-parleur (30 x , 30 N ) et un microphone 

(32 1 , 32 N ) associes, N etant un entier superieur ou 6gal 

a 1, comprenant : 

- des moyens (34) de calcul de signaux auxiliaires, 

30 comportant : 

des moyens (48) de formation de blocs conte- 
nant chacun un nombre predetermine d ' echantillons du 
signal sonore, deux blocs successifs se recouvrant 
mutuellement sur un nombre predetermine d'6chantil- 
35 Ions ; 
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des moyens (50) de ponderation des echantil- 
Ions de chaque bloc par une fenetre d 1 apodisation, 
disposes en sortie des moyens de formation de 
blocs ; 

des moyens (52) de calcul de transformee de 
Fourier, disposes en sortie des moyens de pondera- 
tion ; 

des moyens ( 54 ) de calcul des valeurs de 
frequences inferieures et superieures de bandes 
critiques adjacentes, ou une table de consultation 
contenant ces valeurs ; 

des moyens (56) pour calculer une fonction de 
repartition de 1 1 energie sur les bandes critiques, 
la valeur, dans une bande critique donnee, de la 
fonction de repartition de 1 f energie etant calculee 
en additionnant 1 ' energie. des differentes raies 
spectrales appartenant a cette bande critique ; 

une seconde table de consultation (58), 
contenant des valeurs representant des fonctions 
d'etalement basilaire associees chacune a une bande 
critique donnee ; 

des moyens ( 60 ) pour calculer un spectre etale 
E if en calculant, pour une bande critique donnee, le 
produit de convolution de la fonction de 
repartition de 1' energie et de la fonction d'etale- 
ment associees a cette bande critique ; 

des moyens ( 64 ) de calcul - d 1 un f acteur de 
correction Oi associe a une bande critique donnee ; 

des moyens (66) de calcul d'un seuil de 
masquage frequentiel T ± associe a une bande 
critique donnee, a partir du spectre etale et du 
facteur de correction Oi associes a cette bande 
critique, le seuil de masquage frequentiel T t etant 
defini comme suit : 

10 log 10 T x = 10 log 10 E L - O a - C x 
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ou C A est un parametre supplementaire de correction 
de seuil ; 

des moyens ( 68 ) de generation de signaux 
synthetiques ayant les caracteristiques spectrales 
5 d'un bruit blanc ; 

des moyens (36) de controle, connectes aux N 
voies de signal sonore, pour calculer pour chaque 
bande critique ledit parametre supplementaire de 
correction de seuil C L a partir de plusieurs carac- 
10 teristiques des signaux sonores des N voies ; 

des moyens (70) pour amener le niveau d'un 
signal synthetique a la valeur du seuil de masquage 
frequentiel associe, de fagon a obtenir un signal 
auxiliaire ; 

15 des moyens (72) de calcul de transformee de 

Fourier inverse ; 

des moyens de memoire pour stocker le resultat 
du traitement des blocs successifs ; 

des moyens (74) pour combiner le resultat du 
20 traitement d'un bloc avec le resultat du traitement 

du bloc precedent ; 

des moyens (76) de lecture sequentielle relies 
auxdits moyens de memoire, pour faire passer la 
cadence des donnees de la cadence bloc a la cadence 
25 echantillon ; 

N premiers moyens additionneurs (40 X/ 
40 N ) respectivement disposes sur les N voies de 
signal sonore, pour additionner sur chaque voie le 
signal sonore regu et le signal auxiliaire associe ; 
- N x N moyens de filtrage adaptatif (42 1#1/ . • . , 42 l H , 
•••/ 42 N1/ 42 NN ), comprenant N moyens de filtrage 

adaptatif associes a chacune des N voies de microphone , et 
ayant chacun des coefficients qui forment une estimation de 
la reponse impulsionnelle du chemin de couplage acoustique 
entre le haut-parleur et le microphone associes a une des N 
voies de signal sonore ; 
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- des seconds moyens additionneurs (44 x , . 44 N ), 
connectes en sortie des N moyens de filtrage adaptatif 
associes a chaque voie de microphone, pour additionner les 
signaux de sortie de ces N moyens de filtrage adaptatif ; 

5 - N moyens soustracteurs (46 x , 46 N ), disposes 

chacun en sortie desdits seconds moyens additionneurs (44 x , 
44 N ) et connectes respectivement aux microphones des N 
voies de signal sonore, pour soustraire sur chaque voie de 
microphone le signal obtenu en sortie desdits seconds moyens 
10 additionneurs (44 lf 44 N ) du signal regu par le micro- 

phone de cette voie ; 

- des moyens pour calculer sur chaque voie de 
microphone une erreur d f estimation a partir du resultat 
fourni par le moyen soustracteur connecte au microphone de 

15 cette voie ; 

- des moyens pour corriger de fagon iterative les 
coefficients respectifs des N moyens de filtrage adaptatif 
associes a chacune des N voies de microphone, en fonction de 
1' erreur d' estimation associee a chaque voie de microphone. 

20 6. Dispositif selon la revendication 5, caracterise 

en ce que les caracteristiques servant a calculer ledit 
parametre supplement aire de correction de seuil C ± compren- 
nent le niveau respectif des signaux des N voies et un 
ensemble de valeurs de la fonction d f intercorrelation d'au 

25 moins certaines paires de ces signaux, ledit parametre 
supplementaire de correction C± etant une. fonction crois- 
sante des valeurs de cette fonction d 1 intercorrelation. 

7. Dispositif selon la revendication 5, caracterise 
en ce que les caracteristiques servant a calculer ledit 

30 parametre supplementaire de correction de seuil C x compren- 

nent la somme des energies contenues dans les diverses 
bandes critiques de chaque voie et le rapport des energies 
par bande critique d f au moins certaines paires des differen- 
tes voies, ledit parametre supplementaire de correction C t 

35 etant d'autant plus important que la somme des energies 
contenues dans les diverses bandes critiques est importante 
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et que le rapport des energies par bande critique d'une 
paire de voies est peu different entre les differentes 
bandes critiques, 

8. Dispositif selon la revendication 5, caracterise 
5 en ce qu'il comprend en outre des moyens (62) de calcul d'un 

indice de tonalite a L a partir des energies des raies 
spectrales appartenant aux differentes bandes critiques, 
lesdits moyens ( 64 ) de calcul du f acteur de correction O x 
calculant ledit facteur de correction O t a partir dudit 
10 indice de tonalite a ± . 

9. Dispositif selon la revendication 8, caracterise 
en ce que lesdits moyens ( 64 ) de calcul du facteur de 
correction 0 ± fournissent un facteur de correction 0 i defini 
par : 

15 Oi = max (SO, a L . (k 2 + B ± ) + (1 - a ± ).k 2 ) 

ou max (a, b) designe la plus grande des valeurs a et b, 

ou SO, k 2 et k 2 sont des parametres de valeur predeterminee 

en dB, 

ou a ± est 1' indice de tonalite associe a la bande critique 
20 consideree, et 

ou B ± designe la frequence de la bande critique en Bark. 

10. Dispositif selon la revendication 9, caracterise 
en ce que SO = 24,5 dB, k x - 14,5 dB et k 2 - 5,5 dB. 

11. Dispositif selon la revendication 8, 9 ou 10, 
25 caracterise en ce que les moyens ( 62 ) de calcul de 1 1 indice 

de tonalite fournissent un indice de tonalite constant, 
identique pour toutes les bandes critiques et defini par : 

a t = a = min(SFM/SFM roax , 1) 
ou min(a, b) designe la plus petite des valeurs a et b, 
30 ou SFM pax est un parametre de valeur predeterminee en dB 
associe a un signal sinusoidal pur, et 

ou SFM = 10 log 10 G/A, ou log 10 designe le logarithme decimal, 
G designe la moyenne geometrique de 1 r energie sur un nombre 
predetermine de points de la transformee de Fourier et A 
35 designe sa moyenne arithmetique sur le meme nombre de 
points . 
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12. Dispositif selon la revendication 11, caracte- 
rise en ce que SFM,,^ = - 60 dB. 
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